Possibility of High Tc Superconductivity in doped Graphene 
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Superconductivity at room temperatures, 
with its multifarious technological possibili- 
ties, is a phenomenon that is yet to be realized 
in a material system. The search for such sys- 
tems has lead to the discovery of "high tem- 
perature" superconducting materials such as 
the fascinating cuprates,-'^i2, MgB2^ and most 
recently, Fe-based pnictidesj^ Graphene is a 
remarkable two dimensional conductor f^iSiil A 
newly discovered method to cleave and iso- 
late single or finite number of atomic layers 
of graphene, its mechanical robustness and 
novel electrical properties has caught the at- 
tention of the scientific and nanotechnology 
communities. Undoped graphene is a semi- 
metal and does not superconduct at low tem- 
peratures. However, on "doping optimally" if 
graphene supports high Tc superconductivity 
it will make graphene even more valuable from 
basic science and technology points of view. 
Here we build o n a seven year old sugges- 
tion of iBaskaranl ^ (GB) of an electron corre- 
lation based mechanism of superconductivity 
for graphite like systems and demonstrate the- 
oretically the possibility of high temperature 
superconductivity in doped graphene. 

GB^ suggested the possibility of high tempera- 
ture superconductivity in graphene based on an ef- 
fective phenomenological Hamiltonian that combined 
band theory and Pauling's idea of resonating valance 
bonds (RVB). The model predicted a vanishing Tc 
for undoped graphene, consistent with experiments. 
However, for doped graphene superconducting esti- 
mates of Tc 's were embarrassingly high. Very re- 
cently Black-Schaffer and Doniach^ used GB's ef- 
fective Hamiltonian and studied graphitic systems 
and found that a superconducting state with d -\- 
id symmetry to be the lowest energy state in a 
mean field theory. The mean field theory also pre- 
dicts a rather high value of the optimal Tc- Other 
authors have studied possibiUty of superconductiv- 
ity based on electron-electron and electron-phonon 
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While there is an encourag- 
ing signal for high Tc superconductivity in the phe- 
nomenological GB model, it is important to establish 
this possibility by the study of a more basic and real- 
istic model. Since the motivation for GB model arose 
from a repulsive Hubbard model, here we directly 
analyse this more basic repulsive Hubbard model that 



describes low energy properties of graphene. We con- 
struct variational wavefunctions motivated by RVB 
physics, and perform extensive Monte Carlo study in- 
corporating crucial correlation effects. This approach 
which has proved to be especially successful in under- 
standing the ground state of cuprates, clearly points 
to a superconducting ground state in doped graphene. 
Further support is obtained from a slave rotor anal- 
ysis which also includes correlation effects. Our esti- 
mate of the Kosterlitz-Thouless superconducting Tc 
is of the order of room temperatures, and we also 
discuss experimental observability of our prediction 
of high temperature superconductivity in graphene. 

Low energy electrical and magnetic properties of 
graphene are well described^ by a tight binding Hub- 
bard model defined on a honeycomb lattice with a 
single 2p2 orbital per carbon atom: 
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Here, i labels atomic sites, Cia- is an annihilation op- 
erator for an electron with spin a at site i, Ui^ is 
the number operator at site i of a spin electrons, t 
K, 2.5 eV is the hopping matrix element and U ~ 
6 eV is the onsite Hubbard repulsion. The unique 
band structure of the above model leads to a 'Dirac 
cone' type of spectrum for electron motion close to 
two points in the Brillouin zone, giving rise to a den- 
sity of states that varies linearly with energy near 
zero energy (half filling). 

The Hubbard U is about half the pn free band- 
width, and this places graphene in an intermediate or 
weak coupling regime. Based on this one is tempted 
to conclude that electron correlations are not impor- 
tant. Nonetheless electron correlations are known 
to be important in finite pir bonded planar molecu- 
lar systems such as benzene, naphthane, anthracene, 
caronene etc, all having nearly the same value for 
quantum chemical parameters t and U One of the 
consequences of this is that the first excited spin-0 
state lies above the first excited spin-1 state by more 
than 1 eV. There is a predictable consequence of this 
large singlet-triplet splitting: graphene can be viewed 
as an end member of a sequence of planar pir bonded 
system; this has been also suggested to have a new 
spin-1 collective mode spectrum, as a consequence of 
finite 

Paulingi^ was the first to recognize dominance of 
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singlet correlation between two neighboring pi: elec- 
trons in the ground state. He argued that doubly oc- 
cupied or empty orbitals (polar configurations) 
are less important because of electron-electron repul- 
sion in the 2p2 orbital. Pauling thus ignored polar 
configurations. Once we ignore polar fluctuations 
(states with double occupancy) and consider a reso- 
nance among the nearest neighbor valence bond con- 
figurations we get the well known RVB state. How- 
ever, such a Hilbert space actually describes a Mott 
insulating state rather than a metal. Experimen- 
tally, undopcd graphenc is a broad band conductor, 
albeit with a linearly vanishing density of states at 
the Fermi energy. 

To recover metallicity in Pauling's RVB theory, GB 
combined the broad band feature of pn electrons with 
Pauling's real space singlet (covalcnt) bonding ten- 
dency and suggested^ a low energy phcnomcnological 
model for graphenc: 

where b|^^ = ('^It'^ii ~ ''ii'^it) '^^•^^^'^^ ^ ^Pi^ 
glet on the i—j bond. J (> 0) is a measure of singlet 
or valence bond correlations emphasized by Pauling, 
i. e., a nearest neighbour attraction in the spin sin- 
glet channel. In the present paper we call it as a 
'bond singlet pairing'(BSP) pseudo potential. The 
parameter J was chosen as the singlet triplet split- 
ting in a 2 site Hubbard model with the same t and U, 
J = ((J72-fl6t2)i/2-t/)/2. As [/ becomes larger than 
the bandwidth this psuedo-potential will become the 
famous superexchange characteristic of a Mott insu- 
lator. As shown in^, this model predicts that un- 
doped graphene is a "normal" metal. The linearly 
vanishing density of states at the chemical potential 
engenders a critical strength Jc for the BSP to obtain 
a finite mean field superconducting Tc- The param- 
eter J for graphene was less than the critical value, 
and undoped graphene is not a superconductor de- 
spite Pauling's singlet correlations. Doped graphene 
has a finite density of state at the chemical poten- 
tial and a superconducting ground state is possible. 
B lack- S chaffer and Doniach^ confirmed GB's findings 
in a detailed and systematic mean field theory and 
discovered an important result for the order param- 
eter symmetry. They found that the lowest energy 
mean field solution corresponds to d + id symmetry, 
an unconventional order parameter, rather than the 
extended-s solution. The value of mean field Tc ob- 
tained was an order of magnitude larger than room 
temperature! 

Although results of the mean-field theory are en- 
couraging, it is far from certain that the supercon- 
ducting ground state is stable to quantum fluctua- 
tions. In particular, the GB Hamiltonian does not 
include U which inhibits local number fluctuations. 
In the more basic Hubbard model, a superconducting 
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FIG. 1: Doping dependence of superconducting order pa- 
rameter $ as obtained from VMC calculation of the Hub- 
bard model on a honeycomb lattice for U/t = 2.4. 

state will suffer further quantum mechanical phase 
fluctuations since U inhibits local number fluctua- 
tions. The key question therefore is docs the singlet 
promotion arising out of the local correlation physics 
strong enough to resist the destruction of supercon- 
ductivity arising out of quantum phase fluctuations 
induced by U7 

To investigate the possibility of superconductiv- 
ity in doped graphene we construct a variational 
ground state and optimize it using variational quan- 
tum Monte Carlo (VMC) The ground state we con- 
struct is motivated by the mean-field theory of the 
GB model 

k 

-X(A(fe)(ai^5t_^^-4/_^P + h. c.) 
fe 

where , bf^^ are electron operators on the A and 
B sublattices, fc runs over the Brillouin zone of the 
triangular Bravais lattice, aa,a = 1,2,3 are vectors 
that go from an A site to the three nearest B sublat- 
ticc sites. The free electronic dispersion is determined 
by the function e{k) ~ J2a^^^^" ■ ^^'^ ^^'^ super- 
conducting gap function A(fc) = J2a ^a^^^'^" ■ The 
d+id symmetry motivated by the meanficld solution^ 

i27r(Q-l) 

provides = Ae 3 where A is the "gap pa- 
rameter" , ^/ is a "Hartrce shift" . Starting from this 
mean field theory, we construct a BCS state \BCS)]s[ 
with an appropriate number N of electrons. If we 
work with a lattice with L sites, this corresponds to 
a hole doping of 1 — N/L. Our candidate ground 
state Iv?) is now a state with a Gutzwiller-Jastrow 
factori^iiS g 

\^)=g'^\BCS)N (4) 
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FIG. 2: Dependence of superconducting correlation func- 
tion F{r) on distance r as obtained from VMC calcula- 
tion of the Hubbard model on a honeycomb lattice for 
U/t = 2A,x = 0.2. 



where T) = J2ii''^i-\''^ii + '^iT^'i'i-' operator that 

counts the number of doubly occupied sites. The 
wavefunction Q with partial Gutzwiller projection 
has three variational parameters: the gap parameter 
A, the Hartree shift /i/, and the Gutzwiller- Jastrow 
factor g. The ground state energy {'i'\HH\'^) is cal- 
culated using quantum Monte Carlo methoclii, and 
is optimized with respect to the variational parame- 
ters. Details regarding construction and optimization 
of the wavefunction can be found in the methods sec- 
tion. 

We monitor superconductivity by calculating the 
following correlation function using the optimized 
wavefunctions 



Faj3{Ri — Rj 
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where b^^^ is the electron singlet operator that cre- 
ates a singlet between the A site in the i-th unit cell 
and the B site connected to it by the vector (this 
is just hlj defined earlier with a minor change of no- 
tation). The superconducting order parameter, off- 
diagonal long range order (ODLRO), is 



lim 

\R,-Rj\' 



F{R, -Rj). 
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where F{R, - Rj) = Ea Paa{Rt - Rj)- AH results 
we show in this paper are performed on lattices with 
IS'^ unit cells. 

The superconducting order parameter $ as a func- 
tion of doping, calculated for physical parameters 
corresponding to graphene, obtained using the opti- 
mized wavefunction is shown in fig. [T] Remarkably, a 
"superconducting dome" , reminiscent of cuprates^^ 
is obtained and is consistent with the RVB physics. 
The result indicates that undoped graphene had no 
long range superconducting order consistent with 
physical arguments and mean-field theory^ of the 
phenomenological GB Hamiltonian. Interestingly, 
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FIG. 3: Doping dependence of superconducting order pa- 
rameter as obtained from Slave rotor mean field theory 
of the t-J-U model. 



the present calculation suggests an "optimal doping" 
X of about 0.2 at which the the ODLRO attains a 
maximum. These calculations strongly suggest a su- 
perconducting ground state in doped graphene. 

We now further investigate the system near opti- 
mal doping in order to estimate Tc- Fig. [5] shows a 
plot of the order parameter function F{r) as function 
of the separation r. The function has oscillations up 
to about six to seven lattice spacings and then attains 
a nearly constant value. From an exponential fit one 
can infer that the coherence length ^ of the super- 
conductor is about six to seven lattice spacings. A 
crude estimate of an upper bound of transition tem- 
perature can then be obtained by using results from 
weak coupling BCS theory, using khTc ~ TTei^^" 



Conservative estimates give us kbTc = i.e., is 
about twice room temperature. Evidently, this is an 
upper bound, and an order of magnitude lower than 
the mean-field theory estimates of Black-Schaffer and 
Doniacb^. Further improvement of our estimate of 
Tc becomes technically difficult. It is interesting to 
compare these results with those obtained in a Hub- 
bard model on a square lattice that captures cuprate 
physics. In this latter similar estimate of 

the coherence length ^ is about two to three lat- 
tice spacings2£; however, the hopping scale is nearly 
a magnitude lower and the estimate of Tc is about 
2X1^00111. Again, this provides further support for the 
possibility of high temperature superconductivity in 
graphene. 

It is interesting to see if the results of the VMC can 
be reproduced in a "simple mean-field type" anal- 
ysis that accounts for the presence of U. The re- 
cently developed slave rotor— technique is useful to 
study an intermediate coupling regime that allows 
for number fluctuations at a site, and has been used 
successfully, in the context of cuprates, to study the 
t-J-U model<2^ We adopt a similar approach and use 
a t-J-U model which is equivalent to introducing a 
Hubbard repulsion U in the GB Hamiltonian. The 
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details of the slave rotor technique is presented in 
the methods section. The ODLRO calculated from 
the slave rotor analysis is shown in Fig. [31 and bears 
a remarkable qualitative resemblance with the VMC 
result, i. e., there is an optimal doping that produces 
superconductivity. This result, again, supports the 
possibility of high temperature superconductivity in 
doped graphene. 

It is important to ask about the possibility of 
competing orders that could overshadow supercon- 
ductivity at optimal doping that we have found. 
iHonerkampP ^ has addressed this issue by means of 
a functional renormalization group study of a gen- 
eral Hamiltonian on the honeycomb lattice. He finds 
that in the regime of physical parameter correspond- 
ing to graphene, the system appears to flow towards 
a d + id superconducting state as the temperature is 
lowered. 

It is interesting to contrast the superconductivity 
in cuprates with that in graphene. In the case of 
cuprates (doped Mott insulators), Bloch electrons in 
the entire Brillouin zone are affected by the Hubbard 
U which is larger than the bandwidth. To this extent 
all electrons participate in superexchange or singlet 
bond formation. Consequently correlation hole de- 
velopment is complete, i. e., an electron at a given 
site with an up spin manages to avoid an electron 
with down spin on its site completely, at low energy 
scales. Whereas, in the broad band graphene, only 
those Bloch electrons in the range of energy scale of 
U around chemical potential are affected by the on 
site repulsion. Since this scale is about half the band- 
width, about one half of the electrons are involved in 
singlet bond formation in the ground and low energy 
states. Consequently correlation hole development is 
not complete. The key point is that there is a suffi- 
ciently enhanced singlet correlation, compared to free 
Fermi gas, to be able to support superconductivity in- 
duced by the on site Coulomb repulsions in optimally 
doped graphene. This heuristic picture is supported 
by the variational Monte Carlo result and the slave 
rotor analysis that we have presented in this paper. 

Our prediction of high temperature superconduc- 
tivity raises some obvious questions. Intercalated 
graphite can be viewed as a set of doped graphene 
layers that have a strong 3 dimensional electronic 
coupling. Maximum Tc obtained in these systems is 
around 16 Systems such as CaCg has a doping 

close to optimal doping that we have calculated. Why 
is the Tc so low? On the other hand, superconduct- 
ing signals with a Tc around 60 K and higher have 
been reported in the past in pyrolitic graphite con- 
taining sulfuri^i21 A closer inspection reveals that 
for systems like CaCg i) an enhanced 3 dimension- 
ality arising through the intercalant orbitals makes 
the effect of Hubbard U less important (effect of U 
for a given bandwidth progressively becomes impor- 
tant as we go down in dimensions) and ii) encour- 
agement of charge density wave order arising from 



the intercalant order. Sulfur doped graphite, how- 
ever, gives a hope that there is a possibility of high 
temperature superconductivity. Our present theoret- 
ical prediction should encourage experimentalists to 
study graphite from superconductivity point of view 
systematically, along the line pioneered by Kopele- 
vich and coUaboratorsi^ In the past there has been 
claims (unfortunately not reproducible) of Josephson 
like signals in graphite and carbon based materials 
Again, our result should encourage revival of studies 
along these lines. 

Simple doping of a freely hanging graphene layer 
by gate control to the desired optimal doping of 10 
- 20 % is not experimentally feasible at the present 
moment. It will be interesting to discover experimen- 
tal methods that will allow us to attain these higher 
doping values. A simple estimate shows that a large 
cohesion energy arising from the strong a bond that 
stabilizes the honeycomb structure will maintain the 
structural integrity of graphene. 

The discovery of time reversal symmetry braking 
d + id order— for the superconducting state, within 
our RVB mechanism is very interesting. This uncon- 
ventional order parameter has its own signatures in 
several physical properties: i) spontaneous currents 
in domain walls, ii) chiral domain wall states iii) un- 
usual vortex structure and iv) large magnetic fields 
arising from the d = 2 angular momentum of the 
cooper pairs, which could be detected fj,SR measure- 
ments. Suggestions for experimental determination 
of such an order by means of Andreev conductance 
spectra have been made by iJiang et all^ 

There are also several theoretical and experimen- 
tal issues that needs to be addressed. It is known 
that graphene realized in experimental systems con- 
tains, adsorbed species, inhomogeneities, curvature, 
ripples etci^ Is the superconducting ground state sta- 
ble to these "perturbations"? In particular our the- 
ory gives a substantial ODLRO even for small dop- 
ing. If disorder effects are indeed suppressing a frag- 
ile Kosterlitz-Thouless order in the currently avail- 
able doping regime in real systems, one could uncover 
the hidden superconductivity by disorder control or 
study of cooper pair fluctuation effects. Further anal- 
ysis is necessary to address these issues. 

I. METHODS 
A. VMC 

Diagonalizing the kinetic energy part of eqn. ([3]) 
results in two bands ct and dl . The supcrconduct- 
ing pairing term splits into intra-band pairing and 
inter-band pairing. The latter being unimportant at 
zero temperature can be dropped giving 

Hmf = J2{^^(''K.''k.+E-{k)dljk.) 
k 
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where, Ad{k) = Aq cos (fc • - (p{k)), (p{k) = 
arge(fe). Thus, the variational ground state can be 

written as \BCS) = Uk.^=J^k + "^k'^LZ-kJ^^^ 



where e 



„t 



'fc+cr fccr' k 

ence factors = —a 



'k 



and the BCS cohcr- 



(iJ|)^ + |Arf(fc)P 



To deal with system containing fixed number of 
particles, the variational wavefunction is projected 
onto a fixed number subspace - \BCS)n- The trial 
wavefunction is obtained by applying the Jastrow fac- 
tor to this state - g^\BCS)N- Energy of this state 
is calculated using Monte-Carlo sampling, which is 
then optimized using the simplex method. 



B. SlaveRotor 



The t-J-U model can be written as H 



E 



U J2^ n^^n^i + J S^■SJ. The 



idea is to decompose electronic degree of freedom into 
spinon and charge (rotor) degrees of freedom (DOFs) 
- c]^ = //g.e~'^' where fj^ creates a spin a at site i 
and e***' are the ladder operators for electron den- 
sity, n. 



In this representation, t-J-U model takes the form 
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JJ2{ij) '^i ■ '^j- The total number of parti- 
cles at any site has to be unity to remove unphysical 
states in the expanded Hilbert space. At mean field 
level, the spinon (rotor) DOF can be integrated out 
to give a Hamiltonian in the rotor (spinon) space. 
The two Hamiltonians are coupled via the kinetic en- 
ergy term. They are solved self-consistently using 
standard techniques'^. 
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